Unpredictability of wave function's evolution in 
nonintegrable quantum systems 



in 



> 

ON 



LB. Ivanov 

Theory Department, Petersburg Nuclear Physics Institute, 
188300 Gatchina, Russia 

"<N 
O 
O 

Abstract. It is shown that evolution of wave functions in nonintegrable quantum systems is 
unpredictable for a long time T because of rapid growth of number of elementary computational 
operations 0(T) ~ T a . On the other hand, the evolution of wave functions in integrable systems 
can be predicted by the fast algorithms 0(T) ~ (log2T)^ for logarithmically short time and thus 
there is an algorithmic "compressibility" of their dynamics. The difference between integrable 
and nonintegrable systems in our approach looks identically for classical and quantum systems. 
Therefore the minimal number of bit operations 0{T) needed to predict a state of system for time 
interval T can be used as universal sign of chaos. 
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Chaos as universal phenomenon exists in various systems (for example, in human society) 
and from this point of view the general approach to chaos should not be based on particular 
properties of a system. It is well known |lj that motion of nonintegrable classical systems 
has all attributes of chaos: complexity, unpredictability and randomness. However in nonin- 
tegrable quantum systems the apparent signs of chaos seems to be absent and the main 
direction of studies in the field of "quantum chaos" is semiclassical analysis of various quan- 
tum " signatures" of classical chaos 0, ffl . Nevertheless, taking into account the fundamental 
correspondence principle it is reasonble to suppose that evolution of nonintegrable quantum 
systems should be also unpredictable. As a basis tool to analyze predictability of dynamics 
we use a number of elementary computational operations 0(T) needed to determine a state 
of the system for time interval T. 

Let us consider firstly how a number of elementary computational operations needed for 
prediction of system's evolution depends on time T and accuracy A in classical mechanics. 
It is well known U that a distance between initially close phase space points grows 

with time as 

\\6x(t)\\ = \\8x(0)\\f(t), (1) 

where f(t) ~ e xt for nonintegrable systems and f(t) ~ t for integrable ones. The exponential 
growth of inevitable computational errors leads to unpredictability of long time evolution in 
nonintegrable classical systems. To predict a state (trajectory) of the system by the moment 
T with accuracy A we must make computations with accuracy at least 
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and length of mantissa 

n ~ -log 2 5 = log 2 f(T) - log 2 A. 

It is reasonable to suppose that an algorithm which generates a trajectory has a power 
dependence of number of elementary bit operations 0(n) on the length of mantissa (for 
example, 0(n) ~ n for addition and 0(ri) ~ n 2 for multiplication): 

O(r) ~ n a ~ (Zo<fe/(T) - log 2 A) a , (2) 

where a > 1 is some number. Inserting now f(t) in (^|) we obtain for nonintegrable systems 
(a > 1) 

0(T) ~ T a , (3) 
while for integrable ones 

0(T) ~ (Zo^r)". (4) 

In the first case a number of operations and hence a computational time needed for pre- 
diction grows as a power of time T. Therefore the prediction for a long time is impossible and 
evolution is not algorithmically " compressible" . In regular systems there is a completely dif- 
ferent situation because computational time grows logaritmically with time T and prediction 
is always possible ("compressibility" of evolution). It is interesting to note that algorithmic 
" uncompressibility" of some numerical sequence is the main sign of its randomness: such 
sequence we can only observe. 

Let us consider now a time evolution of some state ip(q, t) in a stationary quantum system 
H(q,p) 

iW t iP(q,t) = HiP(q,t), (5) 
i/)(q,0) = ip (q), 
Any solution of (p3) at t = T can be represented as 

^(g,T) = ^c,e~^ T /^(g), (6) 

where {E^} and {</> M } - are exact eigenenergies and eigenfunctions of Hamiltonian H and 
c ti — ((finl^o)- So if quantum chaos exists it hides in eq. @ which is essence of quantum 
unitary evolution. In analogy with classical case it would be reasonable to investigate a time 
dependence of two initially close wave functions ip\ and %p 2 . It is natural to characterize a 
difference between wave functions by the quantity 



= 11^2-^11 = ^2(l-Re{^ 2 )). (7) 

If at t = we have ||#0(O)|| < e then at t — T the difference (0) is also small ||<5?/>(T)|| < e 
because from unitarity it follows that (ipi\ip 2 ) does not vary with time. Therefore in quantum 
systems there is no divergence of initially close states which is the main sign of chaos in 
classical systems. 

How then unpredictability of evolution can arise in quantum systems? To predict a state 
of the system ip(q,T) at moment t = T it is necessary to know all eigenvalues and eigen- 
functions in (H) exactly. Really we use in (§) approximate eigenvalues and eigenfunctions 
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and therefore by some moment exact wave function and our prediction will be completely 
different. It is obvious that to increase a maximal time of accurate prediction we must 
improve accuracy of eigenvalues. The spectrum of integrable systems is easy computable 
and computational resources grow slowly with increasing of accuracy while for nonintegrable 
systems considerable efforts are needed to improve the spectrum and thus to increase the 
time of correct prediction. Therefore we can assume that quantum integrable systems differ 
significantly from nonintegrable ones in a sense of number of bit operations 0(T) needed 
to predict a future. In the remaining part of this Letter we consider some details of the 
suggested mechanism of unpredictability and hence of chaos in quantum unitary evolution. 

Let us assume that we have found approximate eigenenergies {E^} and eigenfunctions 
{(f) fj}. Then we can predict by means @ that at moment T wave function will be 

^(q,T) = J2~c»e- l ^ T/h Mq), (8) 
n 

where = (4>^\ipo)- Our main goal now is to find a time dependence of the difference 
between exact ip(q,T) and approximate ip(q,T) wave functions 

II^COH = MT)-^(T)\\ = ^2(l-iMMTMT)>)). (9) 
Making use of (@) and (H) we have: 

mnHT)) = j2 c ^ e ~^~ E " )T/h ^M- ( 10 ) 

Now we need some assumptions about accuracy of spectrum obtained. If a maximal error 
of eigenfunctions <^ does not exceed e 

5<frn = 4> tJ ,- <frv, 

ll<%|| < e. 
then the error of is 

5E^ = Ep-Ep = 2E li Re((<f> ll \8<f> li )) + (<%|#|<%), (11) 

\5E^\ < 2\E^\e + e 2 \\H\l 
and error of c M is 

Sc^ = c,j, - c M = (S^^ipo), 

\5Cfj,\ < 6. 

Scalar product {4>n\4>v) we can write as 

(<f>n\4>u) = <V + (</vl<^) = <V + 

\Rp.v\ < £■ 
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Figure 1: A typical behaviour of (a) (ip(T)\ip(T)) and (b) ||#0(T)||. 



Now fllPf) becomes 



Y, c *M 1+R 



iSE, 



" T/h + E c;~c u R^e- l ^- E - +SE ^ h . (12) 



It is clear that first sum carries the main contribution in (i/j(T)\i/j(T)) because second sum has 
a small multiplier R^ u . To understand the main features of ( jID[) it is sufficient to investigate 
time dependence of 



MT)$(T)>«5> M |V 



iSE^T/h 



(13) 



At T = a value of (^(0)|^(0)) equals 1, then (^(^^(T)) decreases and at some moment 
T p reaches values nearly zero. After this (ip(T)\ip(T)) starts to oscillate about zero with some 
small amplitude. For times T > T p the value of ||#0(T)|| oscillates about mean value y/2 
which corresponds to completely independent states. This typical behaviour of (il)(T)\ip{T)) 
and II^COH is shown in Fig. |I|. This means that for T > T p wave function ip(T) completely 
differs from exact function ip(T) and our further prediction is impossible. 

We want now to find how T p and average amplitude of small oscillations A depend on 
errors SE^ and coefficients c M . Our intuition and numerical experience say that T p depends 
on dispersion of errors dE while A is a function of number dim of the most important 
addendums in (|13"D. It is obvious that (ip(T)\ip(T)) becomes nearly zero when arguments of 
exponents in ( |T3"D fill interval [0, 27r]. From this condition we have T p : 

271 nH (14) 



2dE/h dE' 



To estimate A let us make further simplification that all c M in (|T3|) are equal 1/v dim. In 



such a case (fl3T) becomes 



dim 

P(T) = ite(MT)$CT)» = ^— ^cos^T/ft) 



(15) 



and P 2 (T) is: 



P\T) = ^ C os(5E^T/h) cos{5E v T/h) 
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^ £ cos 2 (5E.T/H) + ^~ 2 Y1 cosiSE.T/h) cos(5E u T/h). 



dim 2 ^— ' <iim 2 

It is easy to see that if we make time averaging of P 2 (T) then the first sum equals dim/ 2 while 
the second sum is zero because arguments of cosines are almost independent for different 
and E v . In the result average amplitude of oscillations of Re((ip{T)\ip(T))) can be estimated 

as 



A ~ V p2 ( T ) = "S^=- (16) 



We tested our theoretical estimates (|14]) and (|1^) for wide range of parameters and found 
good agreement with numerical computations. So we have obtained a quite reasonable result 
that uncertainty in our knowledge of spectrum E ± dE leads to limitation in prediction time 
of if>(T) evolution T p ~ h/dE. 

The next step consists of estimating of number of operations needed to obtain the spec- 
trum with some accuracy dE. If the spectrum of a system can be computed by some formulae 
or by effective algorithms (as in the case of integrable systems), the number of operations 
grows slowly with decreasing of error dE. To explain this let us assume for simplicity that 
to compute the spectrum we must perform some series of successive actions which do not 
depend on accuracy required (for example E N = friu(N + 1/2)). Obviously the necessary 
number of bit operations G(dE) depends on mantissa length n ~ —log2dE only 

0{dE) ~ n~< ~ {-log 2 dEy, (17) 

where 7 > 1 - some coefficient characterizing the algorithm of computations. So we have 
obtained that in integrable systems a number of bit operations depends logarithmically on 
accuracy dE. 

Fortunately the variety of systems in nature is not exhausted by easy computable and 
predictable systems. Description of any nonseparable system with strong coupling and es- 
pecially of its excited states is nontrivial problem even for two degrees of freedom. The only 
reliable way to calculate such systems is to use various variation methods, i.e. to minimize 
the energy functional over some space of trial functions. Let us consider how the number of 
operations depends on spectrum accuracy for variational method of Ritz. It is well known 
H that Ritz's method has a power convergence, i.e. approximate eigenf unctions 0^ tend to 
exact ones M as 

\\h-<f>r\\~D- a , (18) 

where D is dimension of space of trial functions and a ~ 1 - some coefficient. Obviously 
(see flllD) the rate of convergence of eigenvalues does not exceed the power law: 



~D~ a . (19) 

To make use of Ritz's method it is necessary to calculate firstly a Hamiltonian matrix with 
number of elements D 2 , and then to find its eigenvalues. The total number of operations to 
do this can be estimated as O(D) ~ D 13 , (3 > 2. To achieve an accuracy e the basis dimension 
D ~ e~ l l a is required (see (0)) and hence the total number of operations is 0(e) ~ E~Pl a . 
Keeping in mind the relation between time of prediction and accuracy of eigenvalues fln|) 
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we obtain the number of bit operations which are needed to predict the wave function's 
evolution over time interval T: 



0(T) ~ T p/a . (20) 
For regular systems (|17|) the number of operations is equal to 

0{T) ~ {log 2 Ty. (21) 



It should be noted that our estimate (20) is valid for any algorithm with power convergence. 

So we see that difference between integrable and nonintegrable systems in our approach 
looks identically in classical and quantum systems. The evolution of integrable systems can 
be predicted by fast algorithms for logarithmically short time (" compressibility" of evolution) 
while in nonintegrable systems such "compressibility" is absent. The reason of unpredictabil- 
ity is universal in classical and quantum systems — rapid growth of number of elementary 
computational operations needed for prediction, but the mechanism of chaos is completely 
different. 
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